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Abstract 



Using the Green function approach to the problem of quantization of the phe- 
nomenological Maxwell theory, the propagation of quantized radiation through 
dispersive and absorptive multilayer dielectric plates is studied. Input-output 
relations are derived, with special emphasis on the determination of the quantum 
noise generators associated with the absorption of radiation inside the dielec- 
tric matter. The input-output relations are used to express arbitrary correlation 
functions of the outgoing field in terms of correlation functions of the incoming 
field and those of the noise generators. To illustrate the theory, the effect of a 
single-slab plate on the mean photon-number densities in the frequency domain 
is discussed in more detail. 



1 Introduction 



It is well known that the use of instruments in optical experiments needs careful 
examination with regard to their action on the quantum statistics of the light 
under study. In particular the presence of passive instruments that may be re- 
garded as macroscopic dielectric bodies responding linearly to radiation can be 
included in the theory by quantization of the phenomenological Maxwell the- 
ory for radiation in linear inhomogeneous dielectrics. The formalism was first 
developed for dispersionless and lossless dielectrics ^ and successfully 

applied to the study of the action of various devices, such as dielectric plates and 
interfaces 0, ^ and optical cavities |]. When broad-band radiation, such as 
an ultrashort light pulse, propagates through (long-distance) dielectric devices 
the effects of dispersion and absorption, which are related to each other by the 
Kramers-Kronig relations, must necessarily be taken into account. In this con- 
text, a number of questions have been arisen, such as the question about the 
velocity at which a single-photon wavepacket travels in dielectric matter under 
the influence of normal and/or abnormal dispersion in the vicinity of an absorp- 
tion line 0, Needless to say that absorption also introduces additional noise, 
at least the (quantum) vacuum noise. 

The problem of describing the effects of dispersive and absorptive linear dielec- 
tric devices on quantized radiation has been considered in a number of papers and 
various quantization schemes have been proposed [T^, |TB|, |T^, |T3|, |TB|, 0, 0. 
As we have recently shown, quantization of the radiation field within the frame 
of the phenomenological Maxwell theory (with given complex permittivity in the 
frequency domain) can be performed using a Green-function expansion of the op- 
erator of the vector potential [[l^, This quantum-field-theoretical formalism, 
which may be regarded as a generalization of the familiar concepts of mode ex- 
pansion, applies to both homogeneous and inhomogeneous dielectric matter and 
is consistent with both the Kramers-Kronig relations and the canonical (equal- 
time) field commutation relations in QED. 

In the present paper we use the method in order to study the behaviour of 
quantized radiation in the presence of a multilayer dielectric plate and to de- 
rive input-output relations. The calculation of input-output relations is com- 
monly based on a development of the familiar formalism of quantum noise theory 
(see, e.g. pO[)' which the explicit nature of the input from a heat bath, and 
the output into it, is taken into account. The advantage of the quantum- field- 
theoretical approach is that it enables one to obtain the input -output relations 
including their dependence on frequency through geometry and dispersive and 
absorptive properties of the layers, because it consistently accounts for the ef- 
fects of radiation-field propagation according to the (phenomenological) quantum 
Maxwell equations. 

It is well known that when dielectric matter in free space is considered and the 
imaginary part of the permittivity may be disregarded (provided that the losses 
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in the chosen frequency interval are sufficiently small), then the input-output 
relations correspond to unitary transformations between input- and output-mode 



operators pl| , p2| . These concepts fail when the effect of absorption of radiation 
through the dielectric matter is taken into account. In this case the outgoing fields 
are not only related to the incoming fields but also to appropriately chosen noise 
sources [|^. Neglection of these supplementary contributions would unavoidably 
lead to a violation of the canonical commutation relations for the outgoing fields, 
and hence effects of quantum noise would left out. Knowing the input-output- 
relations, the quantum statistics of the outgoing fields can be obtained from 
the statistics of the incoming fields and internal (temperature-dependent) noise 
properties of the dielectric matter under study. 

The input-output relations derived may be used in various quantum-optical 
applications, because the multilayer dielectric structure under study may serve as 
a model for a number of devices, such as mirrors, beam splitters, interferometers 
of Fabry-Perot type, and optical fibres. In particular, multilayer dielectric mir- 
rors may be regarded as tunnel barriers for photons and can be used in order 
to measure barrier traversal times 0. The apparently superluminal behaviour of 
photons at dielectric tunnel barriers, which has been observed in recent experi- 



ments, may be modified by the absorption in the dielectric structures These 
and related studies have been a subject of increasing interest in order to answer 
a fundamental question of quantum physics, namely, of how much time it takes 
a particle to tunnel through a barrier. 

The paper is organized as follows. In Sec. the Green function approach to the 
quantization of radiation in dispersive and lossy linear dielectrics is summarized. 
Applying this quantization scheme to radiation propagating through multilayer 
dielectric plates, in Sec. ^ input-output relations are derived. In Sec. | these 
relations are used to study normally ordered correlation functions of the outgoing 
fields. To illustrate the theory, results for the photon-number density statistics 
observed in the case of a single-slab dielectric plate are given. Finally, a summary 
and some conclusions are given in Sec. H. 



2 Quantization scheme 
2.1 Green's function approach 

Let us consider linearly polarized radiation propagating in x direction and allow 
for the the presence of a multilayer dielectric plate characterized in terms of a 
frequency-dependent permittivity e{x,u!) that varies in x direction and obeys, for 
causality reasons, the Kramers-Kronig relations. Introducing the vector potential 
A{x, t) and using the relation 



D{x,t) = eo 



E{x,t)+ f dTx{t-T)E{x,T) 



(1) 
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{E = —A), the classical phenomenological Maxwell equations yield 



A{x,t)-- 



which in the frequency domain reads as 



A{x,uj) = 0. 



dx^ 



(2) 



(3) 



Equation may be transferred to quantum theory as follows [|T7|, p!9|] : 

A{x,uj) =j{x,uj), (4) 



92 ^ ■ 



where the "current" operator 



uo / h 



J[X,UJ) = — 



e,{x,uj) f{x,uj) 



(5) 



is introduced to take into account the noise owing to absorption. In Eq. 
f{x,u) and P{x,uj) are bosonic field operators, 



f{x, u), P{x', u') = 6{x - x') 6{uj - uj'), 



f{x,u;)Jix',cu')] = \p{x,u;)J\x',Lu') 



0, 



(6) 
(7) 



and 



e{x, uj) = e^{x,uj) + i ei{x,uj) (8) 
{A, normalization area in the yz plane). The solution of Eq. (ffl) may be given by 



A{x,uj) = J dx' G{x,x',uj)j{x',uj), (9) 
where G{x, x', u) is the classical Green function obeying the equation 

G{x,x',uj) = S{x - x') (10) 



and tending to zero as x— >-±cx3. It can be shown |T^, |T9|, that the quantization 
scheme outlined above ensures that the (Schrodinger) operators of the vector 
potential and the electric-field strength, 



A{x) 



duj uj) + H.c, 

dujujA{x, uj) + H.c, 



E{x) = 

satisfy the well-known canonical commutation relation 

ih 



A{x),E{x') 



-6{x — x). 



(11) 
(12) 

(13) 
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2.2 Quantum Langevin equations 

From Eqs. (|ll|), ®, and ( p!oD we see that the problem of determining the op- 
erator of the vector potential reduces to the calculation of the classical Green 
function. Before performing the calculations for multilayer dielectric plates, let 
us summarize some results for homogeneous dielectrics, where 



— 1' c/fce^^(^-^') 

271 



2%—n\uj) 
c 



exp 



i—n\uj) |x — X I 
. c 



(14) 



Here 



n\uj) = ^e\uj)= (5{uj)\%^{uj) (15) 

is the complex refractive index. Using Eqs. (^ and (0), Eq. (p!TD may be rewrit- 
ten as 

/?(u; 



A{x) 



dtu 



n 



+H.C., 

(16) 



47rca;eo/5(co')^ n{uj 

where the operators 

a^ix.uj) = -,/27(c^)-e^^(")""/'= dx' e-'"(")""'/77±a;', cu) (17) 

Z V C J-oo 

associated with the waves propagating to the right (+) and left (— ) are intro- 
duced. In the limiting case when the absorption may be disregarded ['y{u!)uj\x — 
x'l/c— s-0], Eq. (|16D reduces to the familiar mode-expansion result fl^, |19]. In 
particular, the operators a±{x, u) and a±(a;, u) become independent on x and sat- 
isfy the well-known commutation relations for photon destruction and creation 
operators, respectively. 

Equation (p!7| ) implies that the operators a±{x, u) and d±{x', u) can be related 
to each other as 



a±{x, ijj 1 



a±(a;',cu)e^^(")"("-"')/^ + 



where 



F±(x,^) = ±-^27(^)-e 



7(x,w). (19) 
Hence, the operators a±{x,uj) obey quantum Langevin equations in the space 

d 



domain 19 



^^a±{x,uj) = =F7(u;)-a±(x,t^) + F±{x,uj), 
where the operators F±{x,uj) may be regarded as Langevin noise sources. 



F±{x,uj),fI{x',uj') = 2-f{uj)-6{x-x')6{u-uj') 



UJ 



(20) 



(21) 



F±(x, u), ai{x', uj')] = ±e {±x' T x) 27(cj)- e-^(")"l"'-"l/^5(a; - uj'). (22) 
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3 Input— output relations 



We now turn to the problem of propagation of quantized radiation through mul- 
tilayer dielectric plates. Assuming dielectric slabs, the interfaces being parallel 
to the yz plane (cf. Fig. 1), the permittivity may be given by 



N 



€{x,uj) = ^ Aj(a;)ej(cc;), 



where 



Xj{x) 



1 if < X < Xj, 
otherwise 



(23) 



(24) 



is the characteristic function of the jth slab (xq— >— cxo, xj^f^oo). In particular, 
for > 3 the system represents an (A^ — 2)-slab dielectric plate surrounded by 
dielectric matter whose permittivity on the left and right, respectively, is ei{uj) 
and eAr(c(j). To determine the Green function G{x,x',uj) that satisfies Eq. (p!0|) 
[and vanishes at infinity], we note that G{x,x',uj) can be decomposed into two 
parts, 

G{x, x', u) = Go{x, x', u) + Gi{x, x', u), (25) 



where, according to Eq. ([T^) , 



N 



Go{x,x',uj) = Xj{x)Xj 



2i—nj{uj] 
c 



-1 



exp 



t — Hjiuj) \X — X 
C 



(26) 



and Gi{x,x',u!) is solution of the homogeneous equation 



dx"^ 



+ ^e(a;,w) 



Gi{x, x' , ijj) = 0, 



(27) 



which implies that 



N 



Gi{x, x', uj) = Y, Aj(x) [C^^\x', ,u)e^"^(^W^ + Cf'ix', ,^)e-^"^('^W^J . (28) 



Clearly, the Cj^^ and Cf {x',uj) must be determined in such a way that 

G{x, x', uj) vanishes at infinity and is continuously differentiable at the surfaces of 
discontinuity. The somewhat lengthy calculations may be performed in a straight- 
forward way. For the simplest case of a single-slab dielectric plate embedded in 
dielectric matter (A^ = 3, cf. Fig. 1), the result is given in [Appendix A] . Because 
of the voluminous formulas, we renounce their presentation for the general case. 

Combining Eqs. and ( P3| ) [together with Eqs. ( PU] ) and (^)], the (Schrodinger) 
operator of the vector potential A{x) for the jth domain (j = 1, . . . , A^) may be 



.(2) 
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represented as, similar to Eq. (pj]), 



roo 



A{x) = du 



^ A7Tcueol3j{uj)A 



nj{uj] 



(29) 

{xj-i<x<Xj), where the dependence on x of the amphtude operators dj±{x, u) is 
governed by quantum Langevin equations of the type given in Eq. (|2y) together 
with Eqs. (|19D and (|2TD [/3{u),'y{u) — > Pj{Lj),'yj{uj)], so that dj±{x,Lj) can be 
represented in the form (0), viz. 



dj±{x,u) = aj-±(x',cj)e^^^(^)"(^-^')/" 

+ rdyF,-±(?/,cu)e^^^(")"("~^)/'= (30) 

J x' 

{xj-i<x,x' <Xj), where 

F,±{x,u) = ±-,/27,(a;)-e^*'3^(-W7'(x,a;). (31) 
t \ c 

Since the relations between the dj±{x,uj) and f{x,uj) sensitively depend on the 
actual expression for G{x,x',uj), the commutation relations (^) [that are based 
on Eq. ([T7|) ] cannot be applied in general. Successive application of Eq. (^) 
enables one to relate the (amplitude) operators di_{x,uj) (— oo <x< Xi) and 
dN^{x,uj) {xn-i<x<oo) forthe outgoing fields to the left and right, respectively, 
to the operators of the corresponding incoming fields, di+{x,uj) and dN-{x,u), 
and operator noise sources associated with the losses owing to absorption. 

3.1 Single-slab dielectric plate 

To illustrate the procedure outlined, let us first study the input-output relations 
for a single-slab dielectric plate of thickness / embedded in two (different) di- 
electrics (A^ = 3, cf. Fig. 1). Substituting in Eq. (|^) for G{x,x',uj) the expression 
(|A 1| ) given in [Appendix A| and introducing the (amplitude) operators di±{x,u!), 
-oo<x<-//2, and a3±(x,cj), //2<x<oo [[Appendix Sll Eqs. ( pi[ ) - ([B^], 



the input operators are found to satisfy the commutation relations 

[ai+(x, cu), dl+{x', Lu')] = e-^^(")l"""'l/^5(cu - cu'), (32) 

[d3-{x,uj),dl_{x',uj')] = e-^3(")l"-"'l/^(5(cu-cj'), (33) 

[di+{x,Lj),dl_{x',u')] = (34) 

[Eqs. ( p 20| ) - ( |B 22| )], so that the input fields from the left and right behave like 
the fields in the corresponding bulk dielectrics and may be regarded as indepen- 
dent variables. Note that di+(x,u;), a3_(x,u;) are commuting quantities. 
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Defining operators 
where 



and 



c±(Z, u) = e 



-72(t^)a;i/c_ 



sinli 



C . 



_|_ g-72('^)wi/c 



and recalling Eqs. 



72 (w) L c J /32(u;) 
) and (|^), we find that 



sm 



c . 



Since 



(35) 
(36) 

(37) 

(38) 
(39) 

(40) 



[Eq. ( |B 23|) ], the incoming- field (amplitude) operators and the operators g±:\uj), (^£-')^(co') 
may be regarded as being independent variables (note that ai+, as- and com- 
mute). Moreover, the g±\uj) and (g^^)"^ {oo), respectively, which are obviously de- 
struction and creation operators of bosonic excitations associated with the plate, 
play the role of the additional operator noise sources in the input-output relations 



(41) 



derived in [Appendix B.l| [Eq. ( |B 5| )], the elements of the 2x2 matrices 



T(i) 



11 ' 
^(1) 



UJ 



12 ' 
(1)/ 



UJ 



ni'iuj) n^'iuj) 



(42) 



and 



Ad) 



(43) 



A^^,\uj) aS¥(^) 

being given in Eqs.( |B 61 ) - ( |B 13| ) [note the simplifications ( [B 14| ) - ( |B 19 ) when 
the plate is surrounded by vacuum]. In Eq. ( p|) the (amplitude) operators of 
the outgoing fields, ai_(— //2,a;) and a3+(Z/2,a;), are expressed in terms of the 
operators of the incoming fields, ai^{—l/2,uj) and d3_(//2,u;), and the operator 
noise sources g^\uj). The characteristic transformation matrix of the plate, T*^^), 
which for an approximately lossless dielectric plate in free space reduces to the 
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well-known characteristic matrix given, e.g., in describes the effects of trans- 
mission and reflection of the input fields, whereas the losses inside the plate give 
rise to an additional matrix, A'-^^ which may be called characteristic absorption 
matrix. 

It should be mentioned that the output amplitude operators di-{x), x< —1/2, 
and 034.(0;), X > 1/2, can easily be obtained from Eq. (0), with x' = —1/2 and 
x' = 1/2, respectively, and application of the input-output relations (Pl). The 
resulting representation of the outgoing fields is of course fully equivalent to the 
Green's function expansion primarily used. The commutation relations for the 
output amplitude operators are given in [Appendix B.l| [Eqs. ( P 24| ) - (p 26|) ]. 
They differ, in general, from those given in Eqs. (|3^) - (|3^) for the input operators. 
The differences vanish when the distances from the plate are large compared with 
the (classical) absorption length or when the plate is in free space. 



3.2 Multilayer dielectric plates 

The results given in Sec. |3.1| can be extended to an arbitrary multi-slab dielectric 



structure (A^>3, cf. Fig. 2) in a straightforward way (for details see [Appendix 
[B.2[ ). In particular, the input-output relations may be given by 



ai-[xi,ijj) 
aN+{xN~i,i^) 



rp(iV-2) 



ai+[xi,uj) 

hN-{xN-l,i^) 



+ A 



(N-2), 

9- (^) 



(44) 



Here, x = x\ and x = xn-i, respectively, are the left and right surfaces of the multi- 
slab plate (note that for the single-slab plate, A^ = 3, the notations X\ = —l/2 and 
xj^-\=X2 = l/2 have been used). The commutation rules for the input operators 
a\j^{x,ijj), ai+(x, cj) (— oo<x<Xi), a^^{x,LS), a}^_{x,LS) (xAr_i < x < 00), and 
the noise operators g^^ ^H^^) are the same as in the preceding section, i.e.. 



\ax+{x,u)),a\^[x' ,uj')\ ■ 
\aj^^{x,uj),a}^_{x! ,J)\ 



-7jv(<^)|a;-i''|/c 



5(cu 



\ax+{x,uj),a}j^^{x' ,J)\ = 0, 



#-^^M,(#~^^(-'))^ =5(.-.'), 



''^iN~2), 

9± 



(.),(#-)(^'))1=o, 



(x,c.),(#-^)(^'))^ 



(45) 

(46) 
(47) 
(48) 
(49) 
(50) 



The input-output relations (0) [together with the commutation relations (^) 
- (pip] apply to arbitrary multi-slab dielectric equipments described in terms of 
a complex permittivity that spatially varies as a multi-step function and whose 
dependence on frequency is consistent with the Kramers-Kronig relations over the 
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whole frequency domain. Typical examples are fractionally transparent dielectric 
mirrors and combinations of them, such as resonator-like cavities bounded by 
dielectric walls. In particular, when the overall device is surrounded by vacuum, 
so that the incoming and outgoing radiation fields propagate in free space, 

ni{uj) = nN{uj) = 1, (51) 

the familiar mode expansions for the input and output (free) fields are recognized. 
For j = 1, Eq. (|2|) takes the form 



A(x) = r du J — - fe'""/'^a.+ (u;) + e-'""/'=a,_(u;)l + H.c, (52) 

^ ^ Jo V 4:TTCujeoA L ■'^^ ' 

where the input and output operators ai±{u!) [al^^uj)] and dN±{u) [a^j^^{uj)], 
respectively, are proper (space-independent) photon destruction [creation] oper- 
ators, and 

The influence of the plate on the incident light through reflection and transmission 
from the two sides are described by the matrix elements ^''(t^), whereas the 
matrix elements A^^~'^\uj) arise from absorption. From Eqs. (^) - (0), the 
bosonic commutation relations for the input-mode operators ai+(cc;), ajsi-iuj) are 
easily seen to be satisfied. Using them and recalling the commutation rules (^) 
- (|50[), the bosonic commutation relations for the output-mode operators ai_(co'), 
aN+{ijj) can then be obtained by means of the input-output relations (p3|) , because 
the relations 




|^(7V-2)|2 . (N-2)^2 I AN-2)^2 , AN-2)^2 

Kll I "T |-'l2 I "T 1^11 I "T 1^12 I 

_ |^(A'-2)|2 |^(A^-2)|2 I 4(A^-2)|2 , AN-2)\2 
~ K21 I "T |-'22 I "T 1^21 I "T 1^22 I 



(54) 



and 



y(7V-2)/y(Af-2)x* ^(Ar-2) .^(7V-2)^ 
11 \ 21 / 12 V 22 



+Asr^'(4r^v + Ar^^(4r^v = o (55) 

are valid [see [Appendix B.l| and [Appendix B.2|| . For notational convenience, the 



frequency argument of the matrix elements T^k'"^"* ^^"^ ^Ifc'"^^ matrices are omit- 
ted. In particular, these relations reflect the fact that the sum of the probabilities 
for reflection, transmission, and absorption of a photon is equal to unity. If the 
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losses inside the plate can approximately be disregarded, A*^^^^) ^ 0, the well- 
known method of unitary transformation is recognized. In this case the relations 
( |5lD and (|55| ) simplify to 

Ti^-'\Ti^-'^r + Tit'\Tit'^r = 0, (57) 

so that T*^^~^) becomes a unitary matrix. Since the photon operators of the 
output and input modes are uniquely related to each other through a unitary 
transformation, the bosonic commutation relations are automatically preserved. 
In general, the T*-^"^) matrix is not unitary and the output-mode operators are 
obtained, according to Eq. (^), from both the input-mode operators and the 
noise operators associated with the losses. The relations and (0), which 
are the natural generalization of the relations ( ^) _ and (|57|) , may be regarded as 
necessary conditions imposed on the T*^^"^^ and A^^"^) matrices of an arbitrary 
dispersive and absorptive multi-slab dielectric device in free space. It should be 
emphasized that these conditions need not be postulated, but they necessarily 
come out of the theory, which also enables one to systematically calculate both 
the f (^-2) and A^^-^) matrices. 



4 Radiation-field correlation functions 



The input-output relations (|53| ) can be used to obtain the quantum statistical 
properties of the outgoing radiation from the properties of the incoming radiation 
and the excitations associated with the dielectric matter. With regard to mea- 
surement, the quantum statistics of radiation is frequently described in terms 
of normally ordered correlation functions, such as correlation functions of the 
electric- field strength or the photon creation and destruction operators (see, e.g., 
[H, |2^). Introducing the notations di = ai_|_, 02 = Oat., d[ = ai_, 0,2 = dj^+ and 
gi=g^^ g2 = g^ from Eqs. (|l^) and (0), the electric-field strength of the 
outgoing radiation in the ith channel (i = 1, 2) reads as 

(58) 




(60) 

{r]i = —1, 772 = 1), where the output-photon operators d[{Lj) can be related to the 
input-photon operators di^u) as, according to Eq. (^31), 

2 



d'Auj) = J2 [Tikiuj) dkiuj) + Aikiu) gk{^^)] (61) 



fc=i 
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[Tik{u>) = T^^ Afc(^) = ^^ik '^\'^))- To express the normally ordered 

electric-field correlation functions of the outgoing radiation, 



m+n 

n E'^i^^^t^ 

/x=m+l 



(62) 



in terms of normally ordered correlation functions of photon creation and destruc- 
tion operators, we use Eqs. (|58| ) - ( |60D and recall the harmonic (exponential) time 
evolution of the photon destruction operators in the Heisenberg picture. We ob- 
tain 



n-m ( ^ \ 



2 



\47rceo^ / 

d^i V^e^'^^^'i ... r du^+n v/^we-'"'"+""'™+-C:';™^f {{u,}) (63) 



{n^='tf^+Vi^x^/c), where 



n 



Co",, 



(64) 



which can be rewritten as, on using the input-output relations (pl|), Eq. 

m 2 



m+n 2 I X 

f_i=m+l fc,i=l 



In particular, when the states of the incoming radiation and the dielectric matter 
are not correlated to each other, the correlation functions of the output pho- 
tons, Eq. (|B3D, can be expressed in terms of sums of products of input-photon 
correlation functions 



k 




k+l 


n ^u^/^) 




n av(^M) 









(66) 



and correlation functions of the excitations associated with the dielectric matter, 



rSS^(K}) 







v+q 






n ^vK) 






ji=p+i 



(67) 



with k,p<m and l,q<n. Clearly, when the matter is prepared in an incoherent 
state, then the correlation functions T}^'^j{{uJ^}) vanish when p ^ q (explicit 
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expressions for the correlation functions observed when the matter is thermally 
excited are given in [Appendix Q) . 

To give a simple example, let us consider the photon-number density (number 
of photons per unit frequency) in the ith output channel (z = l,2), 

K^.i^) = {^H^r<{^)) = c:t'^(^,^) • (68) 

Using Eq. (p3|) and assuming that the overall state is factored, we find that 



k=l 



+ 

where 



T:^{u:)T,2{u:) {a\{u)a2{uj)) + AU^)A,2{uj) {g\{uj)h{^)) + c.c] , (69) 
A^phfc(^) = {al{u)ak{oj)) (70) 

and 

iVdpfc(^) = {gl{uj)gu{to)) (71) 

{k = 1, 2) are the number densities of the input photons and the excitations 
associated with the dielectric matter, respectively. 

Let us first consider the case when the input field is in the vacuum state, 

{a\iuj)a2iu))=0, (72) 

iVphi(a;) = iVph2(^) = 0, (73) 
and the dielectric plate is in thermal equilibrium, 

{glicu)g2ico)) = 0, (74) 

Ndpi{uj)Auj = Ndp2{uj)Auj = - — ntk{uj)Auj, (75) 

zvrc 

where ^ 

"'th(^) = (76) 

M^) - 1 

(T, temperature; k^, Boltzmann constant). For the sake of convenience, in 
Eq. ([75|) we have assumed a large but finite length £ of the quantization vol- 
ume of the radiation, so that the frequency spectrum may be regarded as being 
discrete. Using Eqs. ([72|) - ([75l), we see that Eq. (|^) reduces to 



^ph^(^) = —M^) '^th(^) (t = 1, 2), (77) 



C 

27rc 

where 

= H \ Ak{^^)f (7 



2 



k=l 
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is the (ith-side) absorption coefficient of the plate. The result ( [77D is in full 
agreement with the standard quantum theory of thermal radiation. The plate 
behaves like a thermal radiator which tends to a black body as the absorption 
becomes perfect [ai{u) — >• 1]. Note that for £ — >• oo the photon numbers per 
unit frequency and unit length, N'^^-^^oj) / remain finite. When the plate is in 
thermal equilibrium and the incident light is black-body radiation, 

C 

Npi,i{iu)Aiu = Arph2(cu)Acu = —nthHAuj, (79) 

znc 

the outgoing light is also expected to be black-body radiation. Indeed, combining 
Eqs. (|7l|) , ([7^) , (|7^ , using Eqs. (^), (^), and recalling the relations (|5^ , we 



find that 

Ki^^i^) = ^^^^^i"") (^ = 1,2). (80) 

Finally, let us consider the case when a zero-temperature dielectric plate is 
irradiated from one side, 

iVphi(^)>0, iVph2(^) = 0, (81) 

iVdpi(a;) =iVdp2(^) =0. (82) 

From Eq. (|69D we obtain, on using Eqs. (|7|) and (^ together with Eq. (|8l|) and 



Nl^.iuj) = |7^,i(^)|^iVphi(^) (^ = 1,2). (83) 

The mean photon-number densities in the output channels 1 and 2, respectively, 
are simply given by the mean input photon-number density multiplied by the 
reflection and transmission coefficients. In general, the overall output photon- 
number density is reduced below the input level owing to absorption: 

iV;,,(^)+iV;,2(^)<iVphi(a;). (84) 

In Figs. 3 and 4 the relative photon-number densities of the outgoing ra- 
diation, Afi = N^^,{uj) / Np^iiuj) = |Tn(cu)|2 and A^s = Ar^h2M / ^phi(^) = 
1^21(1^)1^ [Eq. (|83|) ], are shown as functions of frequency and plate thickness for 
a single-slab plate in free space (ei(co') = 63(0;) = 1). The relative photon-number 
density of the radiation absorbed by the plate, A/l = [A^phi(a;) — {Np^i{uj) + 
-^ph2('^))] / ^phi(i^) = «i('^) = «2(<^) is shown in Fig. 5. The results in Figs. 3 - 
5 are given for a simple model permittivity e{uj) = e2{uj) based on a single medium 
resonance, 

eiu) = 1 + , (85) 

which enables one to clearly distinguish the (resonance) region of frequency for 
which the imaginary part of the refractive index may substantially exceed the 
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real part from other (off-resonance) regions for which the imaginary part becomes 
smalL 

For frequencies that are small compared with the medium resonance frequency 
(u;<^co'o) the approximately real refractive index {n{u) = l3{u)+i'y{uj) ^ P{u)>l) 
gives rise to a variation with the plate thickness of the transmitted and reflected 
numbers of photons which show the oscillating behaviour typical for a Fabry- 
Perot device. In this frequency region the losses inside the plate may be dis- 
regarded. Increasing the frequency, both the real and imaginary parts of the 
refractive index, P^uj) and ■yi^uj), respectively, are increased. Increasing •yi^uj) im- 
plies increasing probability for photon absorption in the plate. Since the number 
of absorbed photons increases with the thickness of the plate, the number of 
transmitted photons decreases with increasing plate thickness. Further, owing 
to the increasing absolute value of refractive index the number of photons that 
are reflected are increased at the expense of the number of photons that enter 
the plate. The two effects mentioned become more and more pronounced as u 
approaches ujq. In particular, in the vicinity of the medium resonance the number 
of reflected photons is substantially enhanced. The photons that enter the plate 
are absorbed over a short distance, so that the number of transmitted photons 
rapidly tends to zero when the thickness of the plate is increased. In this "sur- 
face regime" the plate behaves like a lossy mirror, the enhanced reflectivity being 
caused by the large absolute value of the refractive index, which results, in gen- 
eral, from both the real and the imaginary parts. Further increase of frequency 
that is associated with a decrease of the real and imaginary parts of the refrac- 
tive index (region of anomalous dispersion) reduces the effects of strong reflection 
and absorption and the plate again becomes fractionally transparent. Needless 
to say, that for sufficiently high frequencies when the refractive index approaches 
to unity the plate becomes fully transparent. The differences in the response of 
the plate to the incoming radiation for the different regions of frequency are less 
pronounced when the value of the damping constant F in Eq. (|85D is increased. 

In practice, the permittivity is of course a much more complicated function of 
frequency than the model function (^), because of the multi-frequency resonance 
behaviour corresponding to the multi-level excitation spectrum of the medium. 
In Fig. 6 the real and imaginary parts of the refractive index of amorphic Si02 are 
plotted for a region of frequency including absorbing and nonabsorbing sectors 
p7[] . To compare the results for this realistic permittivity with those for the single 
resonance model (|85|) , in Fig. 7 the number of transmitted photons is shown as a 
function of frequency and thickness of the plate. Although Fig. 7 roughly agrees 
with Fig. 4, there are a number of differences in the details which obviously 
result from the multi-frequency resonance structure. We would like to mention 
that, on following the line in Appendix B.2| , the characteristic transformation and 
absorption matrices can be calculated numerically for an arbitrary A^-slab device, 
as we will show in a forthcoming paper. 
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5 Summary and conclusions 



Applying the method of Green-function expansion to the quantization of radia- 
tion propagating through a multilayer dispersive and absorptive dielectric plate, 
we have studied the problem of calculating the proper input-output relations for 
the radiation field operators and presented results for the case when the radi- 
ation propagates perpendicularly to the plate. The plate is described in terms 
of a multistep (spatially varying) complex permittivity in the frequency domain, 
which is introduced phenomenologically and only required to satisfy the Kramers- 
Kronig relations. The advantage of the method is that it enables one to obtain 
input-output relations that do not only apply to regions of frequency far from the 
medium resonances, but they are valid, within the frame of the phenomenological 
linear electrodynamics, in the whole frequency domain. 

In consequence of the inclusion of the losses in the theory the output-radiation 
field operators are found to be related to the input-radiation field operators and 
operator noise sources in the plate associated with the losses, in agreement with 
the dissipation-fluctuation theorem. Disregarding all the losses, the characteristic 
absorption matrix that relates the output-radiation operators to the operators 
noise sources vanishes and the characteristic transformation matrix that relates 
the output-radiation operators to the input operators reduces to a unitary matrix. 
The unitary transformation ensures that the bosonic commutation relations are 
preserved. 

When the multilayer dielectric plate is embedded in an absorbing medium 
the input and output radiation fields can be described in terms of amplitude 
operators whose space dependence (owing to damping) is governed by quantum 
Langevin equations. Only in the limiting case when the surrounding medium 
can be regarded as being lossless (particularly, when the plate is embedded in 
free space), the amplitude operators reduce to the well-known bosonic photon 
operators. In this case, the characteristic transformation and absorption matrices 
of the plate can be shown to be related to each other through conditions that 
ensure preservation of the bosonic commutation relations. These conditions can 
be regarded as the natural generalization of the familiar unitarity conditions for 
lossless plates. 

The input-output relations can advantageously be used in order to obtain 
the quantum statistics of the output radiation from that of the input radiation 
and the noise sources associated with the absorbing matter. With regard to 
measurement, we have considered normally ordered radiation field correlation 
functions. To illustrate the theory, we have studied the input-output relations for 
the (spectral) photon-number densities in more detail. Restricting attention to a 
single-slab dielectric plate, numerical results are presented for a single-resonance 
complex (model) permittivity, and the applicability of the method to realistic 
dielectric matter, such as Si02, is demonstrated. 

Finally, let us briefly comment on the underlying formalism of quantization of 
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the phenomenological Maxwell theory for radiation in dispersive and absorptive 
linear media. It should be mentioned that the formalism used resembles the con- 



cepts of generalized free-field theories [29|. This type of quantum field theory has 
recently been applied successfully to thermo field dynamics for quantum fields 
. The similarities between the two descriptions may provide further insight in 
the basic physical structure of the theory, such as the limit of vanishing absorp- 
tion. In thermo field dynamics this limit corresponds to the zero-temperature 
limit that has been shown to be essentially non-analytical. Similar features are 
also found in an indeterminacy of the dispersion relations. Absorption prevents 
the "spatially damped photons" from exhibiting a well-determined relation be- 
tween energy (frequency) and momentum (wave vector). Similarly, at finite tem- 
perature particle states achieve a continuous spectrum not only for the particle 
momentum but at the same time for the mass parameter (for any fixed value of 
the particle momentum). 
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Appendix A Single-slab Green function 



Let us consider a dielectric slab of thickness I embedded in dielectric matter that 
may be different on the two sides (see Fig. 1). The Green function G{x,x',uj) 
can be obtained from Eqs. (|25D , (^), and (|28D , with = 3. The coefficients 



Cj^\x',(jj) and Cj^'{x',u) in Eq. (pSj) , j = 1,2,3, must be determined from the 
conditions that the Green function is continuously differentiable at the surfaces 
of discontinuity (that is to say, at x = ±1/2) and vanishes at infinity. Note 
that the latter requires the coefficients C[^\x',uj) and {x',uj) to be zero. 
Straightforward but rather lengthy calculation yields [nj=nj{u!)] 



G{x, x') 



2ini — 
c 



e(-x'-i/) {0(-x 



Jin\io\x—x'\l c 



+ 



+ 



2tn2 — 
c 



e (x' + 1/) - e (x' - i/ 
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+ [e(x + i/) -e(a;-i/ 



^gjn.2i^|3::'-(-i)/2|/c^^^gm2a;Z/c _|_ ^in2Uj\l/2-x'\/c^ ^^^^n2U}\l I2~x\l c 
_|_^ |^g«n2a;|x'-(-0/2|/c _|_ ^in2U)\ll2~x' \l c^^^^in-i.u)[ j ^^^^in2Ui\x-{-l) /2\/ c 

+e (a; - \l) d (gi"2c.|x'-(-o/2|/c^^^g 



m2t<j|x' — (— i)/2|/c^ ^in2U>l/c _j_ gm2(^|//2— x'|/c 



^236 



} 



+ 
+ 



2m q 



^ Q(x' - \l) {e (-a; - \l) g»"3c.|x'-V2|/c^g^^gm2C^//c^^^gmi..|(-0/2-x|/c 



6 (a; + i/) - (x - \l) 



^in-j,ui\x' —l/2\/c 



32^9 (e 



,m2aj|//2— a;|/c 



^gin2c.«/c^2ie™2^l^"("')/2|/") + (x - |/ 



^m3a;|a::— x'l/c 



Here, the interface reflection and transmission coefficients Vij =rij{uj) and tij 
Uji^uj), respectively, are defined by 

tij{uj) 



riiioj) — rijiuj) 
rii^uj) + nj{ijj) 

2ni{uj) 



(A 2) 



/ \ , / V ) (A3) 

and the factor = t?(ci;), which arises from multiple reflections inside the slab, 
reads as 



(A 1) 



j=0 



[1 



?^23(t^) 



(A 4) 



Note that the above given form of the Green function permits of a direct phys- 



ical interpretation. The terms in Eq. ( |A 1|) simply correspond to the potential 
propagations of radiation from a source point x' to a point of observation, x. 



Appendix B Photonic operators, noise operators. 

input— output relations 

Appendix B.l Single-slab dielectric plate 



Substituting in Eq. (|^) for G{x,x',u!) the expression ( |A 1| ) and rewriting the 
result (within the space intervals — oo<x<— //2 and l/2<x< oo) in the form 
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(p9[), we easily see that the amphtude operators of the incoming fields from the 
left and right, respectively, ai+(x, cu) and d3_(x,to'), read as 

ai+(x,u) = -./27i(cj)-e-Ti(")""/" T dx' e-'"i("^""'/7"(x', cu), (B 1) 

« V C J-oo 

a3-(x, ^) = \/273(cu)-e^^(")""/" /~ da;' e^^sM-^Vcj^^,'^ ^) (g 2) 
« V c 

[nj(a;) =/3j(u;)+z7j(co')], and the amplitude operators of the outgoing fields to the 
left and right, respectively, ai_{x,u) and ds+{x,u), are given by 



_|_g7i (<^)t<j(a;-//2)/Cg-i?ii {u))uil/c 



ai_(x,cu) = -\/27i(w)-e^i('^)"^'/^ / ' da;' e*"^(")"^''/7(x', cu) 



/?2(^)72(^) 

'^2(c^) \ /?i(t^)7i(u;) 



X 



"^^(u;)+r23(u;)e^"^('^)'^'/^^;(u;) 



+e 



[71 (aj)a::-73(a;)«/2]a;/c 



n3(u;)^ 

^ g-iKH-2n2H+n3H]^V(2c)^^^(^)^^^^^)^(^) a3_(|/,cu), (B 3) 



h+{x,u) = -J273(cu)-e-'^^(-W^ / dx'e-^"«(-)'^^''/7"(a;',cu) 



_^^-^z{u))u){x+l/2)/c^-in-i{u))ujl/c 



r32(^) +t32(^)r2i(^)t23(^)e'*"^('^)"'/^^(^)] a3-(|/,^) 



^ n2(t^) \ /^3(u;)73(a;) 



X 



+e 



- [71 ((^)Z/2+73 ((^)z]aj/c 



^3(0;) 



nl(a;)^J /33(a;) 

^ e-^["i(")-2"2M+«3H]-i/{2c)^^^(^)^^^(^)^(^) ai+(-i/,Lu), (B 4) 

where the operators g'±{uj) are defined in Eq. (|36D . 

Inverting the relations (pS]), we may express the operators in terms of 

the operators g^^\uj) and rewrite Eqs. ( [B 3| ) (for x = —l/2) and ( |B 4[ ) (for x = l/2) 
in the compact form 



^ai_( - ,j,(i)|^ai+( - ( 



a3+(i/,c^) 



a3_(i/,u;) 



(B 5) 
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The elements of the characteristic transformation matrix T^^\ T^^\u), are seen 
to be 



T}l\u;)=e- 



{ui)u!l/c 



X 



ri2(a;) + ti2(cc^)e2"'^--'/V23(c^)^(a;)t2i(a;) 



xi32(c^)e^"^(")"'/=^?(a;)i2iM, 



Til\u;) 



i[/3i(w)+/33(c^)]a;V{2c) 



X 



(B 6) 



(B 7) 



(B 8) 



(B 9) 



and the elements of the characteristic absorption matrix A^^^ A\]^{u}), read as 



1 + e 



72(0^)^'^"^^ g-i/3iHa.V(2c) 



X tuicoMu;)^ fl - e^"2('^)'^'/V23(a;)' 



x72(a;)^e-^3(.M/(2c) 



r2i(uj) + 1 



42^^^) 



72 (C^)"^'^"^^ g-i/33Hc.V(2c) 



X i32(u;)i?(a;)^ [e^"2('^)'^'/V2i(a;) - 1 
Note that when the plate is surrounded by vacuum, 

ni{u) = ^3(0;) = 1, 

so that 

A(a;)=/53(a;) = l, 71(0;) = 73(0;) = 0, 



(B 10) 
(B 11) 

(B 12) 

(B 13) 

(B 14) 

(B 15) 
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the following relations are valid [cf. Eqs. ( |A 2| ) - ( |A 4| )]: 



ri2{uj} =r32{uj) = ^ , ^ = -r2i(u;) = -raal^), 



+ n2{uj) 



321 



t2l{uj) = t21,{uj) 



1 + n2{uj) 
2n2{ijj) 



1 + n2{ijo) 



(B 16) 
(B 17) 
(B 18) 



j=0 



[l -r^iHe 



2in2{u!)wl/c 



(B 19) 



Using Eqs. ( |35D and (^) and recalling the commutation relations (|) and (0), 



from Eqs. ( B 1 ) and (|B 2| ) we derive that 

[ai+(x, w), al+(x', cj')] = e-Ti(^)l^-^'l/'=5(cj - cj'), 

[a3-(x, a;), aUx', u')] = g-^^^-)!^"^' 1/^5 (a; - lo'), 
[ai+{x,uj),al_{x',uj')] = 0, 



a3-(a:',^),(#(^'))^' 



(B 20) 

(B 21) 
(B 22) 
(B 23) 



The commutation relations for the output amplitude operators ai_(x, u), 03+ (x, u) 
and d|_(x, cu), a3_,_(x, cj) can be derived in a similar way or, more easily, using 
Eqs. ( pOD and ( pTD and applying the input-output relations ( |B 5| ) (together with 
the corresponding commutation relations). Straightforward calculation yields 

[ai4x,uj),al_{x',uj')] = 5(cu - cu')|e-^^(")"l"'"'l/^ 



+e 



-ji(cvM\x-(-l/2)\/c+\x'-(-l/2)\/c) 



+i 



7i(^) 

A(^) 



ri(J^(cu)(e'^i(")"' - e-2*/5iH--') - {Til\uj)Y{e 



(B 24) 



[a3+(x,^),at+(a;',a;')]=5(^-^')(e-^^(")"l^-^'l/^ 



+e" 



-rs{i^M\x-l/2\/c+\x'-l/2\/c) 



\Tii\u;)\' + iT2^2^^Hr + lA^VMr + i42^(^)r - 1 



.(1) 



+i 



Tj2V)(e*'^'^'^^''' - e^^^^'^^)'^'^') - (T2^2V))*(e"*^'^'^^''' - e"^*'^^^^)'^^) 



(B 25) 
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[as+ix, iu), al_{x', uj')] = 6(10 - ^')|e-7iH-l-'-(-V2)l/c-73H-|:.-V2|/c) 



X 



^3^^^ |^^-2i/33(a;)c^x7c _ ^-iftHo^i/c 



+ ^(tS^H) -r^ie — -e 



(B 26) 



Needless to say that ai_(x, cj) and a^+{x,uj) are commuting quantities. Using 
Eqs. ( [B 6| ) - ( |B 13| ) and assuming that the plate is embedded in non-absorbing 
media (71 (cj) =73(^7) = 0), the following relations are easily proved correct: 



(B 27) 



T«(.;)(T«H)* + T«(^)(T«H)* 

+aS1)(^)(4V (.;))* + A^^^{u;){AS{co)r = 0. (B 28) 



In this case, both the input-field commutation relations (P 20|) - ( P 22 ) and the 



output-field commutation relations ( |B 241 ) - ( p 26| ) obviously reduce to the famil- 
iar bosonic commutation relations for photon destruction and creation operators. 

Appendix B.2 Multi-slab dielectric plates 

Starting from the results derived in [Appendix B.l| for a single-slab dielectric 
plate (A^ = 3), the results for an arbitrary dielectric plate {N > 3, cf. Fig. 2) 
can be obtained step by step, without explicitly calculating the multi-slab Green 
function. For this purpose we show that when Eqs. - (|50| ) are valid for A^— 1, 
then they are also valid for A^. Using Eq. (Q), with A^ — 1 in place of A^, and 
expressing the operators aN~i±{xN-2, t^) in terms of the operators di±(xi, cu), we 
find that 

( = P f 1 U Q ( €y ] . (B 29) 

\ aN-i~{xN~2,uJ) J \ ai+{xi,uj) J y g'f' -^'(uj) j 

where the matrices P = p(^~3) g^^^^ Q^Q(Af-3) given by 

„ / T^N--i) (N-3)rp(N-3) {N~3)rp{N-3) \ 
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.(7V-3)^(7V-3) JN-3)rj.{N-3) . (JV-3)^(JV-3) . (iV-3)^(iV-3) 

^21 -'12 -'22 ^22 -'12 ~^12 -'22 



-A 



(JV-3) 
11 



-A 



(7V-3) 
12 



(B 31) 



The desired (A^— 2)-slab plate can be obtained by supplement to the (A^— 3)- 
slab plate of the (A^ — l)th slab (cf. Fig. 2). Applying Eq. ( POP to dN~i±{x,Lj) 
{xn-2^x<xn-i) yields 

a;v-i±(x^-i,^) =aiv-i±(x7v-2,c^)e^^^-^("^"'^-^/' + rfVi±(^), (B 32) 
where /jv_i = XAr_i — XAr_2, and 

4_i±(a;) = f dyF^_l±(y,^)e^^--(-)-(^--l-^)/^ (B 33) 

with F^_i±{y,uj) according to Eq. Recalling the commutation relations 

(PID, the operators (ijv-i± found to satisfy the commutation relations 



(iVi±(cj),4-i±(t^')l = '^^{^ - u;')e^^'^"'^"^'^''^"'^'sinh 



7Ar_l(t^)-/Af_l 

(B 34) 



_,^7iv-i(^) 

^ gT«/3jV-l('^)'^(a::jV-2+a;]v-i)/c gjj-^ 



Pn-i[^) 



c 



(B 35) 



By means of linear combination of the operators (i^„^_|_, bosonic operators d]\i_i± 
can be introduced as 



(iAr_i + (u;) 



D 



d'N -i-i^) y 



(B 36) 



dAr_i±(u;),(ijv_i±(t^') = 5(t^-u;')) 

= 0. 



div_i±(u;),d]v_izF(t^') 
In Eq. ( [B 361 ), the matrix D reads as 



D 



11 



21 



|ao| 



— I l«o| 



(B 37) 
(B 38) 



(B 39) 
(B 40) 



22 



Di2 = exp [iarg[ao]] |2 
D22 = exp [2arg[ao]] I 2 



a 



l«o| 



— ) |ao| 



where 



7Ar-.l(u;)-/Ar_l 



sm 



/?iv-i(c^) 

Combining Eqs. (|B 32| ) and (P 36|) we may write 



UJ . 



I3n-i{uj)—In 



where Ru' = Ru 6^, with 



R 



R 



11 



22 ~ ^-^P 



D 



-1 / '^Af-l + (t^) 



-7Af-l(^)-^Af-l 



(B 41) 

(B 42) 

(B 43) 
(B 44) 

(B 45) 
(B 46) 



We now relate the operators d^±{xiy^i,uj) and dN-i±{xN-i,uj) to each other, 
applying Eqs. ( p9D and (|30| ) and recalling that A (a;) is continuously different iable 
at xn-1- Straightforward calculation yields 



/ dN+{xN-l,Uj) 
\ dN_{xN_i,Uj) 

where the elements of the matrix S read as 
^11 ^ 



aAr_l + (xAr_i,t^) 

(In-i-{xn-i,uj) 



\i (3n{uj) 2nN-i{uj) 



S 



12 



-i[/3]V (a;)+/3jv-i(a;)]wx'jv-i/c 



5. 



21 



^ /3Ar(u;) 2niv_i(c^) 



S- 



22 



^ ^Ar(u;) 2raAr_i(u;) 
Combining Eqs. (|B29| ), (|B45| ), and (|B4^ ), we obtain 



(B 47) 

(B 48) 
(B 49) 
(B 50) 
(B51) 



dN+{xN-l,Uj) 

aN_{xN-i,uj) 



SRP 



ai^{xi,uj) 
ai+{xi,uj) 



SRQ 



9 



(iV-3) 



(^) 



SD 



_i / dN-i + {uj) 
dN-i-iuj) 



(B 52) 
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from which we deduce that 



I A' I 9+'^\^) \ ( dN-l + iuj) 



(B 53) 



where the matrices read as 



T(^-2) = {SRP) 



-2 
21 



iSRP) 



22 



iSRP),2{SRP)2i-{SRP)u{SRP)22 (SRP) 



11 



(B 54) 



A' = {SRP)2,' 

^ ( -{SRQ)2i -{SRQ)22 \ 

[ {SRQ)n{SRP)2i-{SRQ)2i{SRP)n {SRQ)^2{SRP)2l~{SRQ)22{SRP)n J ' 

(B 55) 



A" = {SRP)2,' 



21 



22 



V {SD-%,{SRP)2i-{SD-%,{SRP)r, {SD-%2{SRP)2i-{SD~%2{SRP)n J ' 

(B 56) 

Note that g^~^\uj) and d]^_i^{uj) are commuting quantities, because they refer 
to different space intervals [cf. the commutation relations @ or (|2l|) ]. 

Finally, we introduce bosonic operators g^~'^\uj) as linear combinations of 
the operators g'^ ^\^) and dN-i±{^)- 



9+ 

9- M 



u 



, g^+ ^\uj) \ -rJ dN-i + {uj) 



(B 57) 



#~'^M,(^r"^^(^')) 

■#-^^M,(^r'^(-'))^ 



(N-2), 



5{ijj — uj'), 
0. 



(B 58) 
(B 59) 



The elements of the matrix U are given by 



U 



11 



(B 60) 



24 



u. 



21 



U 



12 



22 



where 



exp [iiarg[/io]] |2 
exp [iarg[/io]] |2 

I^'iil' + I^'i2l' + I^nl' + I^' 




I A' 



211 



\A' 



22 



121 ? 
1^21 I ^ 1^221 



4' A'* -^ 4' 4'* -I- A" A"* -^ 4" 4' 

^11^21 ^ ^12^22 ^11^ 21 ^ ^12^ 



22- 



Identifying the the inverse of U with A*^^ 

^ u-1 



(B 61) 
(B 62) 
(B 63) 



(B 64) 
(B 65) 
(B 66) 



(B 67) 



Eq. ( |B 53| ) [together with Eq. ( p 57| )] is the desired result (0). From the struc- 
ture of the Green function [cf. Eqs. (^) - (|28|) ] it is clear that the input operators 
ai+(x,u;) and a7v-(a;, w), respectively, can always be written in the form given in 
Eqs. (|BTD and (p^ [with arbitrary (A^>3) in place of A^ = 3]. With regard to 
the basic-field operators /(x, cu), both the input operators and the noise operators 
g^^ '^\uj) refer to different space intervals. Hence, all the commutation relations 
( ^51) - (|50| ) are satisfied for arbitrary N {N > 3), which implies that the com- 
mutation relations ( [B 24| ) - ( [B 26| ) remain also vahd when d]y^(x,uJ),CL'J^_^_{x',uJ) 

and fW, A«, re- 



and T^N~2)^_^iN-2) ^^^^ substituted for d^^ 
spectively. We finally mention that in a way similar to that outlined above for 
the input-output relations the extension of the relations ( |B 27] ) and ( |B 28| ) to 
arbitrary (A^>3) may be proved correct. 

Appendix C Output correlation functions 

When the dielectric plate is in thermal equilibrium the density operator of the 
matter excitations may be given by 



i=l 



hxjj 



^dp = exp<j - ^ Auj ln[l + nth(^)] 9i{^)9i{^ 



(CI) 



with nth(co') from Eq. (|76|). Note that this density operator of course corresponds 
to the thermal-equilibrium state of the basic field f{x,uj) inside the plate. The 
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correlation functions ( |67[ ) are then calculated to be 



(C 2) 



where the set of indices ( denotes a permutation of the set of indices /i, so that 
uj(^^i <uJt^foT2<(<poTp + 2<(<p + q (note that T}^i'^^{{uJfj,}) vanishes if 
for any ( with p > ( >2 the relation a;^-! = a;^ is fulfilled). 

We now appropriately label the terms that (after disentangling) occur on 
the right-hand side of Eq. (|65D. For this purpose we introduce the set S'^™'"^ of 
arrangements of the m + n indices ( (disposed in ascending order) by assigning 
them to two classes {K = 1,2) and four (possibly empty) groups (j = l,2,3,4). 
The class indices K=l and K = 2 refer to the creation and destruction operators, 
respectively, and the group indices j = 1,2,3,4 are used to distinguish between 
the four excitations to be considered. For = 1 (2) the indices refer for j = 1 to 
d\ (ai), for j = 2 to 4 (^2), for j = 3 to §{ (^1), and for j = 4 to gl (§2). The \f 
indices of the i^th class and jth group are denoted by Cj^(^), i = l, . . . , \f . From 
Eq. (|65[) we then find that 

/(m,7i 

n n n \{%^^^ii^ci{'-))%^,/^e,{P))Ti^2(^^^^^ 

a=l j3=l 7=1 ri=l 1 2 

r) r) 

(C 3) 



c:1Zf ({-.}) 



E c:i:^Sf'^^({a;c.;^cJ)rlt^f({^C3})ri:;f({ 



_5(m,n) 



A} A| A| 



2 ^4, 



n ^vi(.)iHiw) n ^vi,,,2(^ci(x)) 



where 



/^(-^li'^i;-^2'-^2)/ 



'({^CiJ^Ca}) 











-XI 






n 




n 4ki(/3)) 




n^iK?(7)) 




n «2(^ciw) 














_r,=l 



(C 4) 



(the notation Q is used to indicate sets of indices, without distinguishing between 
the two classes). 
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Figure 1: Scheme of the single-slab dielectric plate (2) of thickness / embedded 
in dielectric matter (1 and 3). 
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Figure 2: Scheme of the multilayer dielectric configuration, the arrows together 
with the amplitude operators indicating incoming and outgoing fields. 
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Figure 3: The ratio of plioton-number densities of the reflected outgoing field 
and the incoming field, Afi{uj) = Np^-^^{uj)/Npiii{u), as a function of frequency 
and plate thickness for a single-resonance medium {uq = Ui, T = O.Iuq). 
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Figure 4: The ratio of tlie photon-number densities of the transmitted outgoing 
field and the incoming field, J\f2{uj) = Np^2{^) / ^phi{^) , as a function of frequency 
and plate thickness for a single-resonance medium {uq = ui, T = O.Iuq). 
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Figure 5: The ratio of tlie photon-number densities of the radiation absorbed in 

the plate and the incoming field, A/l = [A^phi(i^) — (-/V'phi(i^)+-^ph2('^))]/-^phi(i^) = 
ai{u!) = a2{uj), as a function of frequency and plate thickness for a single- 
resonance medium {uq = ui, T = O.lcuo). 
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Figure 6: Real and imaginary parts of the refractive index of amorphic Si02 as 
a function of frequency. 
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Figure 7: The ratio of tlie plioton-number densities of the transmitted outgo- 
ing field and the incoming field, N2{oj) = iVpjj2(^)/^phi(t^)) as a function of 
frequency and plate thickness for amorphic Si02 {ujq = 1 ■ 10^^ Hz). 
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